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Abstract. 

The spatial correlation function of the momentum density in the three-dimensional dilute sheared granular fluids is 
theoretically investigated. The existence of the long-range correlation is verified through both analytic calculation and 
numerical simulation. 
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INTRODUCTION 

The time and the spatial correlations play important roles in non-equilibrium statistical physics [1-16]. The behaviors 
of the correlation functions in ordinary fluids are well understood. It is known that there exist long-time tails in the 
time correlation functions in fluids at equilibrium [1, 2, 3,3]. In addition, long-range correlations exist in the spatial 
correlation functions in non-equilibrium ordinary fluids JH 0, . 

On the other hand, interest in the correlations of granular fluids is rapidly growing. In the case of the freely cooling 
state, it is confirmed that long-time tails for the time-correlation of the velocity and the shear stress exist but the 
time-correlation function for the heat flux decays exponentially lloll . In addition, the spatial correlation functions 
are analytically calculated by using fluctuating hydrodynamics and the theoretical results are verified by numerical 



simulations 111 1H . See also the studies on the long-time tails in granular flows 11121 13]. In sheared granular fluids 



lql . Indeed, Kumaran predicted 



there are some studies on time correlations, but the situation is still confusing Il4,ll5l 

that the correlation function satisfies f~ 3c/ / 2 with the spatial dimension d lfl4ll . while we obtained crossover from 
t~ d l 2 to t -( d+7 -)l 2 for the velocity auto-correlation function of nearly elastic granular gases 111511 . However, there is no 
corresponding theoretical or numerical argument on the spatial correlation functions in sheared granular fluids. 

In this paper, thus, we investigate the spatial correlation functions in sheared granular fluids. In section 2.1, we will 
explain the set up. In section 2.2, we will present the analytic results for the spatial correlation of the momentum 
density. In section 2.3, the validity of the analytic results will be tested by our numerical simulation. Finally, we will 
discuss and conclude our results in section 3. 

RESULT 
Set up 

Let us consider a three-dimensional dilute system consisting of N identical smooth and hard spherical particles with 
mass m and diameter a in the volume V. The position and the velocity of the i-th particle at time t are denoted by 
r,-(f ) and v,(f ), respectively. The particles collide instantaneously with each other with a restitution constant e which is 
less than unity for granular particles. When the particle i with velocity v; collides with the particle /' with vy, the post- 
collisional velocities \[ and v' ; - are respectively given by = v ( - — \ ( 1 + e) (n • Vy)n and v^- = Vj + | (1 + e) (n • Vy)n, 
where n is the unit vector parallel to the relative position of the two colliding particles at contact, and v, ; = v, — v ; . 
Let us assume that the uniform shear flow is stable and its velocity profile is given by c«(r) = yy8 a , x , where the Greek 
suffix a denotes the Cartesian component, and y is the shear rate. 

Let us consider the spatial correlation function of the momentum density defined by 

r dr' 

C pp (r)^ — ([p(r + r',0-p(r + r',0c(r + r')]-[p(r',0-p(r' ) 0c(r')]>, (1) 



where p(r,f) = m^f V,-(?)5(r — r/(f)) and p(r,f) = m^f 5(r — r,(f)) are the momentum density and the density, 
respectively. 



Theoretical Analysis 

In order to obtain the analytic expression of C pp (r), we assume that the time evolution of the hydrodynamic fields 
is described by fluctuating hydrodynamics lfl8ll 

d t p+V-(pu) = 0, (2) 
<9,u + u- Vu + V-n/p = 0, (3) 
d,T + u- VT + 2m(n: Vu + V-q)/(3p) = -£r, (4) 

where u(r,f) = p(r,f)/p(r,f) and T(r,t) are the velocity and the temperature, respectively. The heat flux q and the 
pressure tensor Ely consist of two parts as q = q* + q R and = II*,- + IE. Here, q* and IT?- represent systematic 
parts as q* = — kVT — /iVp/m, and ITy — pTSij/m — t] [V/Wj ; + V/W/ — (2Tj/3)5yV^], respectively, where 8ij 
is Kronecker delta. Note that the bulk viscosity disappears in fluids of dilute spherical particles. £, K, ji, and 77 
are the cooling rate, the heat conductivity, the transport coefficient associated with the density gradient, and the 
viscosity, respectively. Here, ji has a finite value when e is less than unity 11711 . q fi and IT^- are the random parts 

of the heat flux and the pressure tensor, respectively. q R and Tlfj are respectively written as qf = Vf^Xf-\ and 
nf; = y/Hjflj, where J* and /?. satisfy (/*) = (/?.) = 0, (^) = 0, (f?(r,t )f s j (r' ) t')) = 25, 7 5(r - r')8(t ~ t'), 

and (ffj( T ,t )n,(r',t')) = 2A ijkl 8(r - r')8(t - 1') with A iJkl = 8 ik 8j, + 8 a 8 jk - 28^/3. 

Let the viscosity r], the cooling rate £, the heat conductivity K", and the transport coefficient associated with the 
density gradient ji be non-dimensionalized as 

77 = 77077*, C = pr£777o, /c=Kb?c*, n = mTK n*/p 1 (5) 

where 77*, £*, JC*, and ^i* are constants which depend only on e in dilute cases [17]. We note 77* ~ 1 and C,* ~ 5e/12 
with e = 1 — e 2 in the limit of small e. Here, 770 and JCo are the viscosity and the heat conductivity in the dilute elastic 
hard-core gas given by 

T7o = aVf, Kb = l5aVf /(4m), (6) 

respectively. The explicit form of the constant a is given by a = 5^Jm/n/{\6o 2 ). 

Here, we introduce the average density po = (p(r,f)), and the average temperature 7b = (T(r,t)). These averages 
satisfy po = mN/V, and 

T = 2m 2 a 2 T 1 *f/(3pfc*)- (7) 

The method to calculate the spatial correlation function is parallel to that for the sheared ordinary fluids @]. We, 
thus, show only the outline of the method in this paper. First, we introduce the fluctuations of the hydrodynamic fields 
8p(r,t) = p(r,f) — po, 8T(r,t) = T(r,t) — T , and 5u(r,f) = u(r,t) — c(r,r). Then, C pp (r) can be approximated by 

f dr' 

C pp {v) = J — ({p + Sp(r + r',t)}Su{r + r',t)-{po + 8p(r',t)}Su{r',t)) 

/dr' f dV 

— <5u(r + r',t) ■ 5u(r',f)) = p 2 J ^C„ M (k) e *' r , (8) 

where we have ignored nonlinear terms of the fluctuations. We have also introduced the correlation function C„ u (k), 
which satisfies 

(5u k (f) • 8u k ,(t)) = (27r) 3 5 3 (k + k')C„ M (k). (9) 
Let us introduce the vector z(r,f ) and its Fourier transform of the fluctuations as 

z r (k,f) = (Sn k (t),ST k (t),Su£\t),8u£\t),Su®(t)), (10) 



where we decompose <5u k (f) as <5u k (?) = 8uj^\t)e^ + du^ (t)t {2) + 8u'£\t)e 1 - 3 ') withe (1) = k/k, e (2) ee {e y - (e (1) • 

e Y )e (1) }/^K, ee e W x e' 2 ', and 8u^(t) = 8v*{f)-e®. Here, e v = (0, 1,0), and jV = \e y - (eW -e y )e^\. Then, the 
linearized time evolution equation for z(k,f ) is obtained from eqs. (|2j, (0, (|4j, and ( TTOb as 

(^-7M^)za(k,r)+ia/j%(k,f) = £ a (k,r). (11) 

Here, the matrix L a| g is given by 

L ap =ikL^+^Lf p +tL^+m ( al+fL ( al, (12) 

where L^l is a matrix depending on £, fc, /x, and rj.R a (k,t) is a random vector which is a function of and IT*. See 

the elements of the matrix and the vector ,R a (k, f ) in Appendix A. 
The solution of eq. (fTTT i is expressed as 

5 ( -t 

z a (k,f)-E/ ds\ j r®(k,t-s)F®(k('Ks-t)),s), (13) 

where 



Y#(k,f) ee ^' ) (k) e -/o^ (i) «W) ) FW(k>0 = ^( k )/? a (M) (14) 
biorthogonal vectors ^ (k), and the eigenvalues A ^ (k) satisfying 



with k(r) ee (fcjcfey — Tk x ,k z ). Here, we have introduced the linearly independent eigenvectors (k), the associated 



+i a/ 0v'i' , (k) = ^ (0 (k)^ ,} (k): (15) 

and^' ) (k)^' ) (k) = ^. 

In order to obtain an analytic expression of \]/$ (k), (p^ (k), and A ® (k), we assume |Ht] 

y-O^/p). (16) 

Here, we note that eq. ( fT6b is rewritten as /~ 2 ~ 0(k 2 ), where Z c is the length scale defined by 



lc=J^. d7) 
V mnoY 



From eqs. (O, (O and 0, the length scale l c becomes 



*3\ 1/4 



^i<r v U^J (18) 

where v is the volume fraction. We, thus, find that l c does not depend on y, but depends on the volume fraction v and 
the restitution constant e because rj* and £* depend only on <?. In the limit of small e, we obtain i c ffv _1 £ - '/ 4 . 
Substituting the solution (fT3l for small A: into eq. (0, we obtain an analytic expression of C pp (r) as 

Cp P (r) =^{Ai(r// c )+A 2 (r/Z c )+A 3 (r// c )}, (19) 



where 



Al (f) = / J**-** r £ /f^(f)[ai-fl 2 cos(c 1 A:a(f))] e - 6 ( r42+,2w '' +;3 ^/ 3 \ 
^ft) 3 Jo 

2%) 



A 3 (f) = / 



(2tt) 3 Jo £ 2 



{2nf Jo k 2 

x [M(k(r)) 2 /fc 2 - 2M(k(t))M(k)Ak(0 +M(k) 2 fe(f) 2 + Jk(f ) 2 ], (20) 



where r = r// c . Here, M(k) = -kk z /{k x ^Jk 2 - k 2 ) tan" 1 {k y /^Jk 2 - k 2 ), a x = (2r/*/3 + fC*/2)/r/*, a 2 = (jc*/2- 

2tj*/3)/tj*, b= (2tj*/3 + k-*/2 + 3^*/4)/tj*, and ci = 2^/5^/^/3^ with v = p 7b/(m?7o(7b)). It should be 
noted that the expressions in eqs. (|T9T > and ( f20b are the same as those for the sheared ordinary fluids by Lutsko and 
Dufty [8]. The difference between ours and theirs exists in the dependence on the restitution constant e through l c . 

Let us explicitly demonstrate the existence of the long-range correlation in C pp (r). Let the angular average of any 
function /(r) be denoted by f(r) = J d£l/(4n)f(r). From the second equation in eq. d20l ), it is easy to show the 
asymptotic behavior 

A 2 (r)°cr- 5 / 3 , r»l. (21) 
Indeed, from the transformation of variables as k' — kr and s — tr~ 2 / 3 , we obtain 

A 2 (r) = [ dk sin ^ f dt ( k2 + 2tk x k y + t2k x) 2 c -(tk 2 + t 2 k t k v +th 2 J3) 
J (2n) 3 kr Jo k 2 

5/3 [Jk^sinW r (r-^k' 2 + 2f- 2 /hk' x k y + s 2 k' 2 ) 2 _ r 4/3 s ^ 2 /3^ + ^ /3) 

J (2tt) 3 k< Jo k' 2 • ( ' 

From this expression, we find that A2(f) satisfies eq. (fJTJ and the integrations in eq. (1221) are reduced to a constant for 
r^> 1 . By using a parallel procedure to that for eq. (f2~TT> . we find 

Ai(r) oc r u/ \ r>l, 

A 3 (f) - f- 5 / 3 , f»l. (23) 
Substituting eqs. (f2TT > and ( f23l into eq. dT9b , we find the long-range correlation in C pp (r) as 

C PP (r) ~ f £ (^)" /3 . r>fe. (24) 
Thus, we confirm that the threshold length l c in eq. (fTTI i plays an important role. 



Numerical Simulation 

To verify our theoretical prediction, we perform simulation of three-dimensional hard spherical particles. In our 
simulation m, a and the average temperature Tq are set to be unity, and all quantities are converted to dimensionless 
forms, where the unit of time scale is awm/%. We adopt the Lees-Edwards boundary condition 11911 . The volume 
fraction v is v c /8 with the closest packing fraction of particle v 6 . We use the parameters y = 0.5 and e = 0.83 to 
keep the temperature unity. We examine the size of the system L = 28,56, 1 12 in our simulation, which contain 4096, 
32768, and 262144 particles, respectively. 

Figure Q] shows the numerical result of C pp (r) for various system size L. There is an apparent finite size effect, 
where C pp {r) decays faster than power-law function for r > 0.3L. For the system with the largest size, we find the 
existence of a region where the spatial correlation function C pp (r) approximately satisfies r~ 5 / 3 . We cannot confirm 
the existence of the region to obey r -5 / 3 in the wide range, but the results seem to be consistent with our theoretical 
result in eq. (l24l . 



CONCLUSION AND DISCUSSION 

Let us discuss our results. From the condition Z^T 2 ~ 0(k 2 ), we could obtain the analytic expression of C pp (r). This 
condition can be rewritten as l c ~ 0(L/(2n)). For the case of our numerical simulation, l c is estimated as 1 .6(7. On the 
other hand, L/ (2n) is estimated as 18(7 in the largest system in our numerical simulation. Hence, it might be suspicious 
that the condition l c ~ 0(L/(2n)) is satisfied in our numerical simulation. However, the result of the simulation seems 
to be consistent with our analytic result. We need more careful consideration on the validity of our analytic method. 




FIGURE 1. The angular average C pp (r) of the momentum density for L = 28,56, 1 12 as a function of the distance r. 



It should be noted that our method in not valid in the true hydrodynamic limit. Indeed, we have assumed eq. ( fT6] l, 
which can be rewritten as ka ~ £ 1//4 . Thus, if e is finite, the applicable range of our analysis is limited for smaller 
systems. We may expect different features in the region of ka <§; £ 1,/4 , which will be our future task. 

The theoretical method we have used is applicable only to the dilute system. In general, sheared granular systems 
in experiments are not dilute ones. Hence, it is not clear that the long-range correlation predicted in this paper is 
experimentally observable. In order to theoretically understand the spatial correlation in such dense sheared granular 
systems, we must improve our theoretical method, which will be also our future work. 

In conclusion, we have analytically calculated the behavior of C pp (r) in three-dimensional dilute sheared granular 
fluids. Based on fluctuating hydrodynamics, we find that there is a long-range correlation in C pp (r). The results are 
verified by numerical simulations. 
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APPENDIX : THE ELEMENTS OF THE MATRIX AND THE VECTOR R a (k,t) 

In this appendix, we explicitly show the elements of the matrix and the vector R a (k,t). Here, we non- 
dimensionalize the variables by the mass m, the characteristic time Vo = po7b/(m77o(7b)), and the characteristic length 
Iq = 2-v/To/w/Vo. The matrix LW and the vector R a (k,t) are explicitly given by 
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k x k y /k 2 

-k x jk L 

kyk^ J kk 


Ri = 


0, 




^2 = 


-^/3jla h ikifi- 


Ri+i — 


■ (') 
a s ie\ ' 








2k x k ± /k 2 

'k-^ky I k 

-k z /k 



(27) 



(28) 



where k\—k 2 ~ k 2 , a s = {mr]* /{pqIq)} 1 ! 2 , an d a h = {5mK* / (3 Po^j)} ■ ff and /?• are the Fourier transforms of ff 
and ffj, respectively. We do not describe the explicit form of and iS 5 ^ in eq. (fT2l because they affect only the 
higher order corrections in eq. ( fT9l ). 
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